Segregation in fluidized versus tapped packs 
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We compare the predictions of two different statistical mechanics approaches, corresponding to 
different physical measurements, proposed to describe binary granular mixtures subjected to some 
external driving (continuous shaking or tap dynamics) . In particular we analytically solve at a mean 
field level the partition function of a simple hard sphere lattice model under gravity and we focus 
on the phenomenon of size segregation. We find that the two approaches lead to similar results and 
seem to coincide in the limit of very low shaking amplitude. However they give different predictions 
of the crossovers from Brazil nut effect to reverse Brazil nut effect with respect to the shaking 
amplitude, which could be detected experimentally. 
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Segregation of hard sphere mixtures is a relevant and 
historically debated problem (see [lj and ref.s therein). 
The phase behavior of these systems is still hotly de- 
bated and it is still controversial whether or not hard 
sphere mixtures segregate in absence of gravity. In the 
last decade a great attention has been devoted to the 
study of the problem of vertically shaken granular mix- 
tures under gravity. It was observed that such systems 
can mix or, under different conditions, segregate their 
components spontaneously according to criteria which 
are still largely unclear, although of deep practical and 
conceptual relevance Rosato et al. |3| demon- 

strated via molecular dynamics simulations that in some 
cases large spheres segregate to the top and small spheres 
to the bottom of the container when subjected to shaking. 
This phenomenon is commonly called Brazil Nut Effect 
(BNE) (while the opposite one, i.e. large spheres on the 
bottom and small ones on the top, is known as Reverse 
Brazil Nut Effect (RBNE) @,0)- The authors suggested 
an interpretation of the BNE as a geometric effect (called 
"percolation" ) where small grains pass through the holes 
created by the larger ones [5, 8]. Along with geometry, 
dynamical effects associated with grains, such as inertia 
|j or convection jjjji which bring large particle up but 
does not allow them to re-enter in downstream, were also 
shown to play a role [Tl). 

Recent results have however outlined that segregation 
processes can involve "global" mechanisms such as "con- 
densation" or phase separation 0, 0|. This sug- 
gested a change of perspective on the issue and the idea 
to formulate a statistical mechanics description of these 
phenomena. In Ref.s 0, [w| a statistical mechanics ap- 
proach was proposed to describe segregation in granular 
mixtures. Assuming that, as in standard thermal sys- 
tems, the grain kinetic energy may play the role of the 
bath temperature, Tt, a thi the system was approximated 
as a standard fluidized gas of elastic hard spheres kept 
at a given temperature, T&at/i, dependent on the shak- 



ing amplitude. This approach may be appropriate to de- 
scribe segregation phenomena in granular mixtures under 
continuous shaking, but not systems subjected to a "tap 
dynamics" (where the energy is pumped into the system 
in pulses, and the measures are performed when the sys- 
tem is at rest). Following the idea, originally suggested 
by Edwards [15| . and recently further studied [la] (see 
[2S for a review), it is possible to develop a statistical 
mechanics description of granular mixtures under such a 
dynamics. Edwards' assumption consists in the hypoth- 
esis that in a granular system under taps, time averages 
coincide with suitable ensemble averages over the "me- 
chanically stable" states, i.e., those where the system is 
found still. Recently, it was shown 0] that this statisti- 
cal mechanics approach holds in good approximation for 
a schematic lattice model of a hard sphere binary mix- 
ture under gravity. In this case two "configurational" 
temperatures have to be introduced in order to describe 
the system's macroscopic states: These temperatures are 
just the inverse thermodynamic parameters canonically 
conjugate to the gravitational energies of the two species, 
and are not related to the grain kinetic energy which is 
always zero. 

In the present paper we study the relation between 
these two different approaches in the framework of a sim- 
ple lattice model for hard sphere binary mixtures under 
gravity. In particular we analytically solve, at the level of 
Bethe approximation, the partition function of the sys- 
tem in two different cases (which we will refer to as case 
I and case II): I) All the configurations are allowed, as 
in a gas (the system is treated as a standard fluidized 
hard sphere gas under gravity, and it is not at rest); II) 
in the framework of Edwards' theory, where the parti- 
cles are required to be in a stable configuration (system 
at rest). We derive the mixing/segregation properties as 
a function of various parameters such as grain masses, 
sizes, numbers and others, and finally we compare the 
results obtained in the two cases. We find that in both 
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cases the system moves from BNE to RBNE increasing 
the mass ratio m%Jm\ or decreasing the concentration 
ratio N 2 /Ni of the two species, instead opposite behav- 
iors are found increasing the shaking amplitudes (i.e. the 
temperatures). In particular in case I the system moves 
from a BNE state to a RBNE increasing the shaking am- 
plitude, and in case II the opposite one is found. More- 
over, for some values of the mass ratio m 2 / mi, the RBNE 
is always found (except at very low temperatures) in the 
case I, while the BNE is always found in the case II. This 
scenario could be experimentally checked. 

The model we consider [l3l Il7j is a hard sphere bi- 
nary mixture made up of two species, 1 (small) and 2 
(large) with grain diameters clq = 1 and V2a Q , under 
gravity on a cubic lattice confined in a rigid box. On 
each site of the lattice we define an occupancy variable, 
nf = 0,1,2, respectively if site i at height z is empty, 
filled by a small or by a large grain. The Hamiltonian is: 
H = Hhc + migHi + m 2 gH 2 , where H x = J2 i>z z8 n -i, 
H 2 = J2i z z $n z 2 are the heights of the two species, and 
TLhc is the hard core potential, preventing two nearest 
neighbor sites to be both occupied if at least one contains 
a large grain. 

The partition function of the system for the case of the 
fluidized hard sphere gas subjected to continuous shaking 
(case I) is given by: 
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where the sum is over all microstates r and ftbath is the 
inverse bath temperature. 

As shown in [l^ the system subjected to tap 
dynamics is well described by Edwards' approach. 
In this case the weight of a given state r is 
0: exp{-n H c(r) - PmngH^r) - (3 2 m 2 gH 2 (r)} II r , 
where (3i and (3 2 are the variables conjugate respectively 
to the gravitational energies of the two species. The op- 
erator IT selects mechanically stable states: II r = 1 if 
r is "stable", else n r = 0. Wc adopt a simple defi- 
nition of "mechanical stability" : a grain is "stable" if 
it has a grain underneath. For a given grain configu- 
ration, r = {rii}, the operator n r has a tractable ex- 
pression: n r = liniR-^oo exjp{-KHEdw} where Hsdw = 

E5„29(5 x-lr.8 z-2„ + <5„z 1 (5 i-l„ I 1 — S z-2„\ 
z,z n i z n. n i n i 1 n. \ n i 2 1 

The system partition function in this case is given by: 

2jj = e -[HHa(r)+0imigH 1 +p 2 magHa] t jj^ ^) 

M 

where the sum is again over all microstates r but, due 
to the projector, only the mechanically stable ones are 
taken into account. We found that the configurational 



temperatures, T\ 
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the limit of low tap amplitudes j3 • Thus for simplicity 
in the following we put l3 con f = Pi = (3 2 . 

Since the exact calculation of Zi and Zn is hardly fea- 
sible, we evaluate the partition functions at a mean field 
level. To this aim we consider a generalization of Bethe- 
Peierls method for anisotropic systems (due to gravity) 
already used in previous papers 0, ^j, i.e., we solve 
the partition function of the system on the Bethe lattice 
shown in Fig by means of recurrence relations (see also 
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FIG. 1: Bethe lattice used in the calculation 

|2*cj). In particular, we consider a 3D lattice box with H 
horizontal layers (i.e., z £ {1, H}) occupied by hard 
spheres. Each layer is a random graph of given connec- 
tivity, k — 1 (we take k = 5). Each site in layer z is also 
connected to its homologous site in z — 1 and z + 1 (the 
total connectivity is thus k + 1). The Hamiltonians are 
the ones above plus two chemical potential terms which 
control the two species concentrations. Hard core repul- 
sion prevents two particles on connected sites to overlap. 

In the following we give the main ideas of the calcu- 
lation in the case I, and refer to a longer paper |l8j in 
preparation for the details, and the calculation in the case 
II (see also [HE3)- The lattice of Fig. His iterated in 
three different directions, "side" , "up" and "down" . As 
consequence three different branches exist, and the corre- 



sponding partition functions, Z^fJ ~ and Z^'^ ud ^ (re- 
stricted respectively to configurations in which the site 
i is empty or filled by a particle of species n = 1, 2) 
can be defined. The Bethe-Peierls recursion equations, 
which allow calculation of the partition function, are 
more easily written in terms of the local "cavity fields" 
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e J " =Z { j£/z(% z) , (with n= 1,2). 

The fluid phase corresponds to a solution of Bethe- 
Peierls equations where local fields in each layer are site 
independent |2*H . Such a solution, characterized by hor- 
izontal translational invariance, is given by the fixed 
points of the following equations: 
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as function of the tap amplitude, and tend to coincide in 
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where n = 1, 2; ^2 are the chemical potentials of the two 

species and mi,a are the grain masses; An^ = l+5 n ie h i ' , 

- 1 + e fc i" + Ui z) = 1 + S nl A'\ D { n z) = 1 + 

6nie f i"\pW = l+e 9 ™ +e 9 ' z) andQ^ = l+e f ^ +e^"\ 
From the local field the system free energy can be derived: 

F = f)AFW-^f; AJrW-^A^W, (4) 
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A similar calculation can be developed in case II. In the 
following instead of using the chemical potential vari- 
ables, \i\ and (j,2, we will use the conjugate variables, 
iVi and N2, the number per unit surface respectively of 
the small and large grains. 

From the free energy F, we calculate the density pro- 
files, 171,2(2), defined by: 
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where N is the site number on each layer, n = 1, 2 corre- 
sponds respectively to small and large grains, and T\ and 
Ti are respectively the radius of small grains (ao/2) and 
of large ones (v / 2ao/2). In Fig. 0we compare 171,2(2) ob- 
tained in the two cases for low and high temperatures. At 
low temperatures the density profiles are quite similar, as 
expected since, in the case I, the system explores almost 
exclusively the mechanically stable states. At high tem- 
peratures the density profiles are instead much different. 

In Fig. AH/H = (i?i - H 2 )/(H 1 +H 2 ) is plotted as 
a function of the concentration, N%, of the small grains, 
at a given concentration of the large ones, N 2 : In both 
cases the system shows a crossover from BNE to RBNE 
as the concentration is decreased. 

In Fig. 0]a diagram of the vertical segregation is given 
in the plane mass ratio - temperature (respectively Tbath 
and T con f in the two cases) for fixed values of Ni and N 2 . 
The broad lines correspond to the crossover from BNE 
to "mixing" , and from "mixing" to RBNE (the "mixing 
region" is defined as -0.1 < AH/H < 0.1). The lines 
on the left are those relative to the case I, which ap- 
proach asymptotically the vertical line m 2 jm\ = 1 in the 
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FIG. 2: (color online). Density profiles in case I, where all the 
configurations are allowed (blue stars and dashed lines), and 
in case II, where only the stable configurations are allowed 
(red empty circles and continuous lines) for mi = 7712 = 1 
Ni — 3.6, N2 = 2.8, and respectively T CO nf = 0.5 and T bat h = 
0.5 (top), and T con / = 5 and Tbath = 5 (bottom). In all the 
cases small grains are on the box bottom (BNE). 




FIG. 3: (color online). AH/H as a function of Ni (for N 2 = 
0.4) in the case I (blue continuous line on the top) for Tbath = 
2, and in the case II (red dashed line on the bottom) for 
T con f = 2. The broad horizontal lines identify the "mixing 
region" where -0.1 < AH/H < 0.1. 



limit Th a th —> 00, as in a point-like bidisperse gas under 
gravity. At high temperatures in both cases the system 
moves from BNE to RBNE with a continuous crossover. 
Increasing the temperatures the mixing regions broaden 
since the entropic term in the free energy becomes more 
important, and the curves in the two cases get apart since 
the differences of the entropic terms increase. At very low 
temperatures in both cases the mixing region becomes 
a line due to the presence of a real phase transition in 
the corresponding model without gravity 13]: small and 
large grains tend to demix and due to the gravity one 
of the two species is energetically favored to stay on the 
bottom depending on the mass ratio. By further decreas- 
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FIG. 4: (color online). Diagram of the vertical segregation 
state in the plane mass ratio - temperature (respectively Tbath 
and T conf ) for ATi = 0.8 and N 2 = 0.4. The blue lines on the 
left are relative to the case I, and the red lines on the right 
are relative to the case II. 

ing the temperatures the curves in the two cases become 
closer, and finally tend to a common value in the limit 
T ^ 0. 

It is evident that at given mass ratio the crossover from 
BNE to RBNE is obtained by varying the shaking ampli- 
tude (i.e. the temperatures) in opposite ways in the two 
cases: In the case I a crossover from BNE to RBNE is 
found when Tb a th is increased as experimentally observed 
in , and in the case II a crossover from RBNE to BNE 
is found when T con f is increased. Moreover, for some val- 
ues of the mass ratio, the RBNE is always found (except 
at very low temperatures) in the case I, instead the BNE 
is always found in the case II. This scenario could be 
checked via molecular dynamics simulations and experi- 
ments of binary granular mixtures: Increasing the shak- 
ing amplitude of the continuous external shaking (case 
I) and of the tap dynamics (case II) in fact corresponds 
respectively to increasing Tb at h and T con f. 

In conclusion we focus on the problem of the vertical 
size segregation in binary granular mixtures subjected to 
external driving. We compare the predictions of two sta- 
tistical mechanics approaches proposed to describe gran- 
ular materials, the first more appropriate to deal with 
fluidized systems subjected to continuous shaking and 
the second with granular materials subjected to tap dy- 
namics. Applying the two approaches to a simple hard 
sphere lattice model under gravity we find, in contrast to 
similar features at low shaking amplitudes, very different 
behaviors (which could be experimentally checked) as the 
shaking amplitude is increased. This is due to the pre- 
vailing of two different mechanisms giving rise to segrega- 
tion: for the fluidized granular mixture the segregation is 
mainly energetic driven, instead for the tapped mixture 
the segregation is mainly entropic driven. A consequence 



is for example that the BN region in the second case is 
larger than in the first one, due to the fact that the stable 
configurations with small grains on the box bottom are 
more than those with large grains on the bottom. 
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